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Abstract 

We investigate combined effects of topology and boundaries on the vacuum expectation 
value (VEV) of the fermionic current in the space with an arbitrary number of toroidally 
compactified dimensions. As a geometry of boundaries we consider two parallel plates on 
which the fermion held obeys bag boundary conditions. Along the compact dimensions, 
periodicity conditions are imposed with arbitrary phases. In addition, the presence of a 
constant gauge held is assumed. The nontrivial topology gives rise to an Aharonov-Bohm 
effect for the fermionic current induced by the gauge held. It is shown that the VEV of 
the charge density vanishes and the current density has nonzero expectation values for the 
components along compact dimensions only. The latter are periodic odd functions of the 
magnetic flux with the period equal to the flux quantum. In the region between the plates, 
the VEV of the fermionic current is decomposed into pure topological, single plate and 
interference parts. For a massless held the single plate part vanishes and the interference 
part is distributed uniformly. The corresponding results are generalized for conformally-flat 
spacetimes. Applications of the general formulas to finite-length carbon nanotubes are given 
within the framework of the Dirac model for quasiparticles in graphene. In the absence of the 
magnetic flux, two sublattices of the honeycomb graphene lattice yield opposite contributions 
and the fermionic current vanishes. A magnetic flux through the cross section of the nanotube 
breaks the symmetry allowing the current to flow along the compact dimension. 

PACS numbers: 03.70.+k, ll.10.Kk, 61.46.Fg 



1 Introduction 

In quantum field theory, the boundary conditions imposed on the field operator give rise to 
a number of interesting physical effects. In the first class of models, these conditions are due 
to presence of boundaries having different physical nature, like macroscopic bodies in QED, 
extended topological defects, horizons and so on. In this type of problems the field operator 
obeys the boundary condition on some space- like surfaces. A well-known quantum effect induced 
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by this kind of boundary conditions is the Casimir effect (for reviews see Ref. pQ). The Casimir 
effect is among the most striking macroscopic manifestations of nontrivial properties of the 
quantum vacuum. The boundary conditions imposed on the field operator alter the zero-point 
modes of a quantized field and shift the vacuum expectation values (VEVs) of quantities, such 
as the energy density and stresses. This leads to the appearance of forces acting on constraining 
boundaries. 

In the second class of models, the boundary conditions arise due to the nontrivial topology of 
the space. The latter induces periodicity conditions imposed on the field operator along compact 
dimensions. This type of models appear in many high-energy theories of fundamental physics, 
including supergravity and superstring theories. Models of a compact universe with nontrivial 
topology may also play an important role by providing proper initial conditions for inflation in 
the early stages of the Universe expansion [2]. An interesting application of the field theoretical 
models with compact dimensions appeared in nanophysics recently. For long-wavelengths, the 
dynamics of quasiparticles in a graphene sheet is well described in terms of the Dirac-like theory 
in 2-dimensional space with the Fermi velocity playing the role of the speed of light (see Ref. [3]). 
In the geometry of a single- walled carbon nanotube, which is generated by rolling up a graphene 
sheet to form a cylinder, the background space for the corresponding Dirac model has topology 
R 1 x S . For toroidal carbon nanotubes, an additional compactification along the tube axis 
leads to the topology (S 1 ) 2 . 

In models with nontrivial topology, the periodicity conditions along compact dimensions give 
rise to Casimir-type contributions in the VEVs of various physical observables. In Kaluza-Klein 
models, the topological Casimir effect has been used for the stabilization of moduli fields and 
as a source for dynamical compactification of extra dimensions. As it has been discussed in 
Refs. [3], the Casimir energy related to the compact subspace of extra dimensions can serve 
as a model for dark energy driving the accelerated expansion of the universe at the present 
epoch. Note that, the recent measurements of the Casimir forces between macroscopic bodies 
provide a sensitive test to constrain the parameters of various types of long-range interactions 
predicted by unification theories [5]. The influence of extra compactified dimensions on the 
Casimir effect in the geometry of two parallel plates has been recently discussed for scalar [6j, 
electromagnetic [TJ, and Dirac fermion [8j |9] fields. The combined quantum vacuum effects 
from boundaries (branes) and from the non-trivial topology of spatial dimensions in braneworld 
models on anti-de Sitter bulk are considered in Ref. [10] (for the Casimir energy and stresses in 
braneworlds sec Ref. |llj). In these models the Casimir forces provide a natural mechanism for 
stabilizing the radion field, as required for a complete solution of the hierarchy problem between 
the gravitational and electroweak mass scales. 

In the papers cited above, the vacuum energy and the forces acting on the constraining 
boundaries were considered. For charged fields an additional important characteristic of the 
vacuum state, bilinear in the field, is the VEV of the current density. In Ref. [12], the VEV of the 
current density is evaluated for a massive fermionic field in spaces with toroidally compactified 
dimensions. It has been shown that the nontrivial topology of the background spacetime leads 
to the Aharonov-Bohm effect for the fermionic current induced by the gauge field. Continuing in 
this line of investigation, in the present paper we consider the combined effects of boundaries and 
compact spatial dimensions on the VEV of the fermionic current. Although the corresponding 
operator is local, due to the global nature of the vacuum, this quantity carries an important 
information about the topology of the background space. The current acts as the source in the 
Maxwell equations and therefore plays an important role in modelling a self-consistent dynamics 
involving the electromagnetic field. As the boundary geometry here we will consider two parallel 
plates on which the fermionic field obeys the MIT bag boundary condition. In addition, we 
assume the presence of a constant gauge field (the fermionic current in a conical space with a 
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circular boundary has been recently discussed in Ref. [IS])- 

The outline of the paper is as follows. In the next section, we specify the mode functions 
and the eigenvalues of the momentum for the Dirac equation in the region between two plates 
with bag boundary conditions on them. The VEV of the fermionic current in this region is 
investigated in Sect. [3j This VEV is decomposed into a pure topological, single plate and 
interference parts. We also give a generalization for conformally-fiat background spacetimes with 
toroidally compact dimensions. In Sect. H] we apply the general formulas for the investigation 
of the fermionic current in finite-length metallic and semiconductor carbon nanotubes, within 
the framework of the Dirac model for quasiparticles in graphene. The main results of the paper 
are summarized in Sect. In Appendix we derive an alternative expression for the VEV of the 
fermionic current in the geometry of a single plate. 

2 Mode functions 

We consider a Dirac fermion field i/j(x) in background of a (D + l)-dimensional flat spacetime, 
in the presence of a gauge field A p . The dynamics of the field is described by the Dirac equation 

{i-fDp - m)V(x) = , (1) 

with D p = d p + ieA^. The dimensions for Dirac matrices 7^ = (7 ,7) are Nd x Nd, where 
Nd = 2^ D+1 ^ 2 ^ and the square brackets denote the integer part of (D + l)/2. In the discussion 
below we use the Dirac representation: 

with cr = (<7i, . . . ,o\d). For the matrices a p in Eq. ([2]) we have the anticommutation relations 
a n a v + a v (X t l = 2<5jm/. In a two-dimensional space the irreducible representation corresponds to 
Nd = 2 and the corresponding Dirac matrices can be taken as 7 M = (<rp3, icpi, icp2), where op M 
are the Pauli matrices. 

We assume that q = D — p — 1 spatial dimensions with the Cartesian coordinates z q = 
(z p+ 2, ■ ■ ■ ,zd) are toroidally compactified. The length of the Z-th compact dimension will be 
denoted by Li, so ^ z\ ^ L\ for I = p + 2, . . . , D. The remaining coordinates z p+ i = 
(zi, . . . , z p+ i = z) are not compactified with —00 < z\ < 00, I = 1, . . . , p + 1. Hence, we consider 
the background space with topology R p+1 x Along the compact dimensions we impose 

on the field operator quasiperiodic boundary conditions: 

^(t,z p+1 ,z q +Liei) = e 2niai ip(t,z p+1 ,z q ), (3) 

with constant phases \ai\ ^ 1/2 and with being the unit vector along the direction of the 
coordinate zi, I = p + 2, . . . ,D. For untwisted and twisted fermionic fields one has special values 
0.1 = and cti = 1/2, respectively. 

Our main interest in this paper is the VEV of the fermionic current = ip'j^vp, where 
-ip = ip^° is the Dirac adjoint and, as usual, the dagger denotes Hermitian conjugation. We 
assume the presence of two parallel plates placed at z = and z = a, on which the field operator 
obeys MIT bag boundary conditions 

(1 + vfn„) i>{x) = , z = 0,a, (4) 

where is the outward oriented normal to the boundary. From these conditions it follows 
that on the boundaries n p j^ = 0, i.e. the component of the current along the normal to the 



3 



boundaries vanishes on the plates. In what follows we will consider the region between the plates, 
< z < a, with = — <5^ +1 and = <5^ +1 for the plates at z = and z = a, respectively. The 
expressions for the VEVs in the regions z < and z > a are obtained by the limiting transition. 

The VEV of the fermionic current is expressed in terms of the two-point function S$ (x,x') = 
(0\[ip r (x),ip s (x')]\0) , where r, s are spinor indices and |0) stands for the vacuum state. The 
expression for the VEV reads: 

{f{x)) = (Otf"(*)|0> = -±Tr(<fsW(x,x)). (5) 

Let {tp^~\x),tp^ \x)} be a complete set of positive- and negative-energy solutions to the Dirac 
equation with a set of quantum numbers f3 specifying the modes (see below). Expanding the 
fermionic field operator in terms of these solutions with the coefficients being the annihilation 
and creation operators, the VEV of the current is presented in the form of the mode sum 

(f) = -E^W^r^) -$ + w^ +) (*)]- (6) 

This expression is divergent, hence some regularization and subsequent renormalization proce- 
dure is necessary. The important point here is that, owing to the flatness of the background 
spacetime, the structure of divergences is the same as for the topologically trivial Minkowski 
spacetime. As a result, the renormalization is reduced to the subtraction from the VEV of the 
corresponding Minkowskian quantity. In what follows, in order to make the expression on the 
right-hand side of Eq. ([6]) finite, we will assume the presence of some cutoff function, without 
writing it explicitly. The special form of the latter will not be important in the remainder of 
the discussion. 

For the evaluation of the VEV by Eq. Q we need to have the mode functions (x) for 
the equation ^Q) obeying the periodicity conditions ([3]) and the boundary conditions (j4|). For 
the gauge field we assume that = const. Although the corresponding field strength vanishes, 
the nontrivial topology of the space leads to the Aharonov-Bohm-like effect for the fermionic 
current. In the problem under consideration, the mode sum ([6]) involves the integration over the 
momentum along uncompactified dimensions and the components of the constant gauge field 
along these dimensions are simply removed by shifting the integration variable. Consequently, 
the VEV does not depend on these components. For that reason we will assume a nonzero 
vector potential along compact dimensions only: = (0, —A) with A = (0, A q ) and A q = 
(A p+ 2, • • • , Ad). 

With this choice, the mode functions are obtained from those discussed in Ref. [9] by a simple 
generalization: 

*P = 4^ ( • ieA) ^ ) , (7) 



For the upper and lower components of the positive- and negative-energy functions one has the 
expressions 

p(±) = gi^-il fr&eikp+i* + ^Wg-'W), ( 8 ) 

with Z|| = (zi, ... , z p , z p+ 2, • • • i z d) being the coordinates parallel to the plates. For the corre- 
sponding momentum we have ky = (k p , kg ) and k p = (fcj, . . . , k p ), = (A;^ 2 , . . . , &^)- 
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Now, the energy corresponding to the modes reads 



(±) = Jk2 + k 2 p+1 + (kf } =F eA,) 2 + m 2 . (9) 



From the periodicity conditions ([3]), for the eigenvalues of the momentum components along the 
compact dimensions we find 

k\ ±] = 27r( ni ± ai )/L h I = p + 2, . . . , D, (10) 

with ni = 0, ±1,±2, .... For the components along the uncompactified dimensions one has 
— oo < ki < oo, I = 1, . . . ,p. 

The relations between the coefficients in Eq. ([8]) is found from the boundary condition @ 
on the plate at z = 0: 

, s k p+1 a p+ ia^ • (ki +) - eA) - m(u; (+) + m) - k^ +1 , +) 

ip\ = -; t— 7 r (f_ , 

{m — ik p+ \) [lo^ + > + m) 
f _l Vfio'p+i "!! • ( k || + eA)-m(a;(-) + m) - k 2 +l 

r- = " — -< \-u v (-) _i_ ^ —v+ > 11 

where cry = (<t 1( a p , a p +2, ■ ■■> cd)- As we have already extracted the coefficients in 

Eq. (JT|), we can impose the normalization conditions tp_ (f_ = 92+ <p+ = 1- Now, from 

Eq. (jlip it can be seen that one has also the relations cp± + (p+ = ^v 7 - ^ = 1- As a set of 

independent spinors with this normalization we take ip^ = and ip+ ^ = w^' . Here w^ a \ 

(7 = 1,..., Nd/2, are one-column matrices having Np/2 rows with the elements u>, = 5i a: and 
w (°)' = iw^. 

It remains for us to impose on the mode functions the boundary condition (|4|) on the plate at 
z = a. From this condition we find the following equation for the eigenvalues of the component 
of the momentum normal to the plates: 

ma sin(k p+ ia) / (k p+ ia) + cos(/c p+ io) = 0. (12) 

We denote the positive solutions of this equation, arranged in the ascending order, by — 
k p +ia, n = 1, 2, — Note that, for a massless field, one lias \n = 7r(n — 1/2). Hence, the mode 
functions ([7j) are specified by the set j3 = (k p , kg^, n, a). 

The coefficients in Eq. © are determined from the normalization condition with the 

integration over the region between the plates. They are given by the expression 



A (±)2 _ ^ (=fc) + m 
p " 4(27r)PLU^aV q 



sin(2A n ) 



2X n 

with V g = L p+ 2 • • • Ld being the volume of the compact subspace. 



(13) 



3 Fermionic current 

Having the complete set of mode functions, we can evaluate the VEV of the fermionic current 
by using the mode sum formula ([6]). First of all, we consider the /x = component which 
corresponds to the vacuum charge density. 
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3.1 Charge density 



By taking into account Eq. (|7|) for the positive- and negative-energy mode functions, for the 
charge density one finds the expression below: 



N D 

4aV a ^ 



ma 



dk 



v 



(2vr)P 



E 



n=l 

e 2i\ n z/a 



1 



sin(2A n ) 
2A„ 



-2i\ n z/a 



+ 



ma — i\ n ma + iX 



(14) 



where n g = (n p+ 2, • • • ,^d), — °° < n i < +oo. This expression is not convenient for the direct 
evaluation of the charge density. This is related to the fact that the eigenvalues A n are given 
implicitly and the terms with large values of n are highly oscillatory. These two problems are 
solved by using the Abel-Plana-type summation formula |14^ [T5] 



71=1 



TTf(Xn) 



sin(2A n )/(2A n ) 



7rma/(0) 
2(ma + 1) 



+ 



dx f(x) 



/(-«) 



x+ma 2x 
x—ma 



+ 1 



(15) 



The corresponding function f{x) in Eq. (|14|) is an even function and, hence, the last integral in 
Eq. (|15p vanishes. In addition, one has /(0) = 0. Further, for the part in the first integral in 
the right-hand side of Eq. (|15|) we have 



dx 



— + 

m + ix 



2ixz p +i 



m 



ix 



0. 



(16) 



Hence, after applying Eq. (|T5|) we get: 



0'°> 



iVp 
2aV„ 



dk 



P+i 



(2vr)? 



(17) 



The latter coincides with the charge density in the spacetime with topology R p+1 x (5 1 ) 9 , in the 
absence of boundaries. The renormalized value of the latter vanishes (12j . Hence, we conclude 
that the presence of boundaries with MIT bag boundary conditions does not induce any vacuum 
charge density. 



3.2 Current density 

Next we consider the spatial components for (j^), \i = 1, ... ,23. Firstly, let us consider the 
component of the current along the direction normal to the plates. From Eqs. (|6|) and (|7|) the 
following expression is obtained 

• _ A U)2 

o ,p+1 ) = JEE -J—frM+^ip® - ( w i)4 >°' } ]- (is) 

3=+- P 

By making use of Eq. ([8]) , it can be seen that 

^ )+ d p+1 ^ ~ (d^upM+frV) = 2ik p+1 (^ )+ V f - = 0. (19) 

Hence, we conclude that the VEV of the normal component of the fermionic current vanishes. 
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Now we turn to the components of the current along the directions parallel to the plates. 
From Eqs. © and © we get 



o" 



E E 



A 



(J)2 



(?) 



m 



W+ Wi^)]. 



where Z = 1, 
notations cr 



(-) 



+ 2, . . . D, j' = + for j = — and f = — for j = +. We have also used the 
= and cr^ = <7+. By using Eqs. (J8|) and (fTTj) . after lengthy calculations, the 
following representation is obtained for the corresponding mode sum in the region between the 
plates: 



4aV„ 



E 



dk 



(27T) 



k, 



X ■ 



n=l 

e 2i\„z/a 



1 



ma 



+ 



sin(2A n ) 
2A n 

e -2i\ n z/a 



-1 



where 



and 



with k„ 



u; \ ma — «A n ma + iA n 

fcj = 2-K(m + a/)/L z for Z = p + 2, . . . , D, 



UJ 



(k 



p+2, 



'A2 /a 2 + k 2 + k 2 + m 2 , 
, Zee). In Eq. (f2T|) we have defined 

Qi = eAiLi/{2ir) - a t . 



(20) 



(21) 



(22) 



(23) 



Hence, the VEV of the fermionic current depends on the phases in the periodicity conditions 
and on the components of the vector potential in the combination (|23p . From Eq. (120|) it 
follows that the VEV is a periodic odd function of this parameter with the period equal to 1. 
In particular, the fermionic current is a periodic function of A[L\ with the period of the flux 
quantum $0 = 2vr/|e| ($0 = 2irhc/\e\ in standard units). 

From Eq. (|20j) it directly follows that the components of the current density along the 
uncompactified dimensions vanish: 



(j 1 ) = for I = 1, 



,P- 



(24) 



We also can see that the fermionic current along the Z-th compact dimension vanishes for the 
special cases = and di\ = 1/2. In particular, this is the case for untwisted and twisted 
fermion fields in the absence of the gauge field. 

By the reasons given after Eq. (|14p , formula (|20p is not convenient for the investigation of the 
properties of the current. In order to obtain a more workable representation, for the summation 
of the series over n we use Eq. (fT5|) . As a result, the VEV of the current is decomposed as: 



(j 1 ) 



-I)(0) + (/)(!) 

(x 2 - q 



2irV n 



q n q GZi 



m 



-1/2 



X 



x+m 2ax 
x—m 



+ 1 



m 



dk p 
(2vr)P 

,2zx 



dx 



-2zx 



+ 



m — x 



m + x 



(25) 
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where k?, = k 2 + k 2 , . Note that the first two terms in the right-hand side of Eq. (|'25p come from 
the first integral in the right-hand side of Eq. (|15p . In Eq. (|25p , 

W « n %J ( 2 ^) P+1 Jk2 +1 + k 2 + m 2 



is the corresponding fermionic current in the boundary-free spacetime with topology R p+l x (S 11 ) 9 
and the part 



(J 



Iv (1) N D m 



/dk \ f°° 



e 2lxz /(m — ix) + e~ 2lxz /(m + ix) 



x 2 + k 2 + k 2 + m 2 



is induced by the presence of a single plate located at z = 0, when the second plate is absent. 
The latter can be seen by the direct evaluation of the VEV for a single plate, or by taking the 
limit a — > oo in Eq. (|25p . Hence, the last term in Eq. ()25[) can be interpreted as being induced 
by the presence of the second plate at z = a. 

The boundary induced parts in Eq. (I25p are finite (see below) and the renormalization is 
necessary for the boundary- free part only. The latter has been investigated in Ref. [12J. The 
corresponding renormalized VEV is given by the expression 



U 



MO) 



2N D m D+1 L l ^ 
(2 ^(D+i)/2 2^n lS m{2^n m ) 



„,_S,-. ^' ' ° 5 - l) (m 9 (L 5 ,n,))<-«)/ 2 ' (28) 

where n g _i = (n p+2 , . . . ,raz_i,n m , . . . ,n D ), a g _i = (a p+ i, . . . ,aj_i, 3j+i, . ..ao), if^z) is the 
Macdonald function (i.e. the modified Bessel function of the second kind) and 



5 (L 9 ,n,) = ^ L^J . (29) 

The reason for the sign difference of ol\ in the expression of 6l\ in Eq. (|23p and in the corresponding 
expression of Ref. |12] is that, in the latter reference, for the evaluation of the VEV the negative- 
energy modes have been used, with the eigenvalues fcj (see Eq. (|10p ) instead of k\ . This 
means that, in fact, the formulas given in Ref. are for the periodicity conditions © with ai 
replaced by —a/. The formula (|28p is obtained by making use of the zeta function technique. An 
alternative representation, derived by using the Abel-Plana-type summation formula, is given 
in Rev. [12]. In the discussion below we will be mainly concerned with the effects induced by 
the plates. 



3.2.1 Fermionic current induced by a single plate 

First of all, we discuss the part in the fermionic current induced by a single plate. As it is seen 
from Eq. (|27p . f° r a massless field the latter vanishes. In the case of a massive field, for the 
further transformation of Eq. (|27l) . we rotate the integration contour in the complex plane x by 
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the angle n/2 for the term with e 2lxz and by the angle —tt/2 for the term with e 2lxz . This gives 
the result: 

, m _ N D ,n E I J- «-*V(m + ,) 



where 



"in, = yk 2 + m 2 . (31) 
This expression is further simplified by using the relation 

V A/i5TS:./^-k2- m i r(( P +i)/2)7 t 



The final formula can be written in a form valid for both plates and for both sides of the plate: 

"°° —2x\z—a,j 



, , A poo c —2x\z—a i \ 

(j l )f ] = ~^N D m Y,h dx (x 2 - ml q )<*-W e , (33) 



with the notation 

(47r) -( P +l)/2 

r((p + i)/2) 

In Eq. ([33]) . j = 1,2, ai = 0, 02 = a, and is the VEV induced by a single plate placed 

at z = a,j, so (j 1 )^ = (j 1 )^ ■ Note that the integral in Eq. f)33[) is a monotonically decreasing 
positive function of m nq . An alternative expression for the plate induced part in the VEV of 
the fermionic current is derived in Appendix [Aj by using the summation formula (I62p . 

If the length of the r-th compactified dimension, L r , r ^ I, is large compared to the other 
length scales, the dominant contribution to the sum over n r in Eq. (|33|) comes from large values 
of n r . In this case, to the leading order, we can replace the corresponding series by an integral, 
using the relation 

+00 j ^oo 

J2 /(27r|n r + a r |/L r )^^ / dyf(y). (35) 
The double integral is reduced to the single one with the help of the formula 

' dy r dx (x 2 -y 2 - 6 2 )(p-D/ 2 /(*) = ^j^^P [°° du(u 2 - b 2 f' 2 f(u), (36) 



2r(l+p/2) 

and from Eq. (|33p the corresponding expression is obtained for the topology i? p+2 x (5 1 ) 9-1 . 
Note that this way of calculating does not work for the case r = I, as the leading term obtained 
by the replacement (|35p vanishes. The asymptotic expression for large values of will be given 
below, by using the alternative representation of the boundary induced part in the VEV given 
in Appendix. 

In the opposite limit, i.e. when L r is small, the behavior of the current density depends 
crucially whether the parameter a r is zero or not. For a r = the dominant contribution in Eq. 
(|33p comes from the mode with n r = 0, and to the leading order one has 

(j i )f ] » a^a^!^)- v>5Vx ( ^-> ( 37 ) 
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where (j }\ RP+ i x , s i\ q -i is the corresponding quantity in (D — l)-dimensional space with topology 
KP +1 x (S 1 ) q ~ 1 and with the lengths of the compact dimensions L p+2 , ■ ■ ■ , L r _i,L r+1 , . . . , Lp. 
For a r 7^ 0, the dominant contribution comes from the region near the lower limit of the 
integration in Eq. (133j) and we get: 

As before, the contribution of the term with n r = dominates and the VEV is exponentially 
suppressed by the factor e ~^ a r\z-a^ |/L r _ ]\[ t e that this asymptotic form is valid for r = / as 
well. 

For the investigation of the boundary induced part near the plate and for large values of L;, 
it is more convenient to use Eq. (|65p . From this expression it is seen that the boundary induced 
part is finite on the plate and the corresponding value is obtained by the direct substitution 
z = 0. This contrasts with the case of the fermionic condensate and the VEV of the energy- 
momentum tensor (see Ref. [9j) which are divergent on the boundary. 

In the limit of large values of L\ the dominant contribution in Eq. (|65p comes from the 
region near the lower limit of the integral over x. Assuming \at[\ < 1/2, we can see that the 
contribution of the term with = 0, i = p + 2, . . . , I — 1, 1 + 1, . . . , D, dominates in the sum 
over n g _i. To the leading order we get 

~ NpLi sin(27ra;) P /2+i m , Ll r , q v 

with the notation 

D 

m\= ^TTOi/Lif + m 2 . (40) 

Hence, for large values of the length of the compact dimension we have an exponential suppres- 
sion. In this limit the total VEV is dominated by the boundary-free part (j 1 )^ ' ■ As we see, 
in both limits of small and large values of L; the boundary induced part in the VEV of the 
fermionic current goes to zero. 

In the numerical examples below we consider the simplest Kaluza-Klein model with a single 
compact dimension with p = D — 2, q = 1, and k nq = 2n\nD + au\/Li). In the left panel of 
Fig- Q3 the current density induced by a single plate at z = is plotted versus the parameter 
old, for a fixed value mz = 0.5. The numbers near the curves are the corresponding values of 
mLo (the length of the compact space in units of the Compton wavelength of the fermionic 
particle). As we already mentioned before, the current density is a periodic function of the 
parameter an with the period equal to 1. As it is seen from the graphs, the absolute value of 
the parameter qd, for which the current is maximum, increases with increasing the length of 
the compact dimension. The value of the current at the maximum decreases with increasing the 
length. The right panel of Fig. [T] presents the dependence of the fermionic current on the length 
of the compact dimension for separate values of the parameter old (numbers near the curves) 
and for mz = 0.5. As we already mentioned, the boundary induced part vanishes in both limits 
of small and large values of the compact dimension length. 



10 




Figure 1 : VEV of the current density induced by a single plate in the Kaluza-Klein model with 
spatial topology R 3 x S 1 , as a function of the parameter (left plot) and of the parameter 
mLi£) (right plot), for a fixed distance from the plate corresponding to mz = 0.5. The numbers 
near the curves are the values of mL^ (left plot) and old (right plot). 



3.2.2 Geometry of two plates 

By using Eq. (|32p . the total VEV in the region between two plates is presented in the form 

= (3 l ) {0) + (J l ) (1) ~ N D ^ £ h dx 



if 



( x 2_ m 2j(p-l)/2 
x+m c 2ax _|_ J 



2 — m 



n 2zx 



-2zx 



+ 



m — x m + x 



(41) 



where m nq is defined by Eq. (|31|) . Combining with Eq. (|33|) for the geometry of a single plate, 
we find 



{j i } (°)-2N D ^ y: h r 

q n„eZ9 Jm »> 



dx 



{x 2_ <q) ( P -m f me 



x+m 2ax i i 

x—m " 



X — m 



cosh [(a — 2z)x 



(42) 



As we could expect, this expression is symmetric with respect to the plane z = a/2. The integral 
in Eq. (|42p is a monotonically decreasing positive function of m nq . 
The formula (142 j) is further simplified for a massless field: 

A r oo ( x 2 _ U2 \(p-l)/2 

o-') = o-') (0) - E k > dx n " 

9 v. e7,<J •'few. 



where 



e 2ax + ! 



(43) 



£ = 



(44) 



i=p+2 



Note that in this case the distribution of the fermionic current in the region between the plates 
is uniform. An alternative expression is obtained by making use of the expansion (e y + 1) _1 = 
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— ^Z^ = i{—^) n e ny ■ Then the integral is expressed in terms of the Macdonald function and one 
gets 



= (J l ) {0) + 



AN D a- p ' 2 
(4vr)P/ 2 + 1 y (! 



E ^< 2 E 



q n„eZi 



n=l 



—^K pl2 (2nak nq ) 



(45) 



At large separations between the plates, compared with the lengths of the compact directions, 
the dominant contribution to Eq. (j45]) comes from the terms with rii = for |dj| < 1/2 and 
from the terms rii = 0, ±1 for on = ±1/2. For |d;| ^ 1/2, to the leading order one finds: 



(46) 



where /?q = X]£ P +2 (^i/^i) 2 an d -^1/2 i s the number of compact dimensions for which on = 
±1/2, i 7^ I. As it is seen, in this limit the boundary induced part in the VEV is exponentially 
small. 

Extracting the parts corresponding to the single plates, the fermionic current can also be 
presented in the form 

= {j i)(o) + -£ {j i)f) +A (j l ), (47) 
where the interference term is given by the expression 



A P N D 
V„ 



n q eZi 

rn 
x + m 



00 ( x - 
dx — 



m 



2 \(p-l)/2 



x+m c 2ax 
x—m 



+ 1 



e ~2zx _|_ e 2x(z-a) 



(48) 



For a massless field the interference part coincides with the boundary induced part. In the limit 
of large separation between the plates, the dominant contribution comes from the region near 
the lower limit of the integration in Eq. (|48|) . To the leading order we get 



2 N WnN D onm<* 1)/2 m - m c _ 2amQ 
K(47ra)(P +1 )/ 2 L; m + m e 



771 



mo + m 



e -2zm , e 2m (z-a) 



where 



D 



(2na i /L i ) 2 + 



m 



(49) 



(50) 



i=p+2 

For a massless field this result is reduced to Eq. (|46p . Hence, in the limit under consideration, 
the interference effects between the boundaries are exponentially suppressed. 

Up to now, we have considered the region between the plates, < z < a. For the regions 
z < and z > a, the VEV of the components of the fermionic current along the compact 
dimension z% is given by the expressions: 

(J 1 ) = (j l ) {0) + (j 1 )?, z<0, 

(J 1 ) = (j^ + ti 1 )?, z>a, (51) 

with the plate induced part given by Eq. (f33j) . In particular, the boundary induced part vanishes 
in these regions, for a massless field. 
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Figure 2: Boundary induced part of the fermionic current versus old (left plot) and L^ja (right 
plot) in the region between two plates for the model with spatial topology i? 3 x 5 1 and for a 
massless field. The numbers near the curves are the values of a/Ln (left plot) and of old (right 
plot). 

As in the case of a single plate, in the numerical example we consider the model with spatial 
topology R 3 x S 1 and with the length of the compact dimension Ljy. In Fig. [21 the boundary 
induced part in the VEV of the corresponding fermionic current is plotted for a massless field in 
the region between two plates, as a function of the parameter old (left plot) and as a function 
of Ld/ a (right plot). The numbers near the curves correspond to the values of the ratio a/Ljj 
for the left plot and to the values of old for the right plot. Similar to the single plate case, 
the boundary induced part vanishes in both limits of small and large values of the compact 
dimension length. 

3.3 Fermionic current in conformally-fiat spacetimes 

From quantum field theory in curved spacetime it is well-known (see, for instance, Ref. |16j ) 
that, in conformally-flat spacetimes, the expectation values of bilinear field combinations for 
conformally invariant fields are related to the corresponding quantities in flat spacetime by a 
conformal transformation. By taking into account that a massless fermionic field is conformally 
invariant in all spatial dimensions, we can use the results given above for the generation of the 
topological and boundary induced parts in the VEV of the fermionic current for conformally-flat 
spacetimes with spatial topology R p+l x {S l ) q and with the line-element 

D 

ds 2 = n 2 ( Zl )[dt 2 -J2(dz t ) 2 }, (52) 
i=i 

where, as before, ^ Z\ ^ Li, I = p + 2, . . . , D. The functions £l(zi) are assumed to be periodic 
along the compact dimensions. The Dirac matrices in this spacetime, 7m) > are related to the 
matrices in flat spacetime by 7^ = e^7 ( , where ef = _1 <5^ is the corresponding tetrad field. 
If iprn) is a massless fermion field in the spacetime with the line-element (|52j) . then we have the 
following conformal transformation: t/jfQ) = Q~ D ^ 2 ip, rfcft) = ^~ D ^ 2 ^P- By taking into account 
that, for the normal to the boundaries, one has JVn)^ = ^n^, we see that the bag boundary 
condition [1 + ^Tm) 71 ^)^]^^) = is invariant under the conformal transformation. Assuming 
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the presence of two boundaries at z = and z = a (z p+ i = z), with bag boundary conditions, 
for the fermionic current in the spacetime (|52p . we get 



(53) 



Here, the first term in the right-hand side is the corresponding VEV for a spacetime with the 
line-element (|52p and with spatial topology R D , and the second term is induced by the nontrivial 
topology and by the boundaries. With this form, the renormalization is needed for the first term 
only. For the special cases of de Sitter and anti-de Sitter (AdS) spacetimes, one has ft 2 = l/(Ht) 2 
and Q 2 = l/(kz) 2 , respectively. 

Let us consider the case of AdS bulk in detail. The boundaries in AdS spacetime can play the 
role of branes in higher-dimensional generalizations of Randall-Sundrum-type braneworlds with 
toroidally compactified internal spaces. In these models the coordinate z is compactified on an 
orbifold S 1 /Z2, and the branes are on the two fixed points. The left/right boundaries correspond 
to the hidden/visible branes. By taking into account that z = represents the AdS boundary 
[z = 00 corresponds to the AdS horizon), we will assume that the boundaries are located at 
z = a\ 7^ and z = 02, a\ < 0,2- The physical coordinate in the direction perpendicular to the 
branes corresponds to y = fc _1 ln(/cz), with k~ l being the AdS radius. The physical distance 
between the branes, in terms of this coordinate, is given by yo = ln(a2/ai). To discuss 
the physics from the point of view of a L>-dimensional observer residing on the visible brane, 
we introduce rescaled coordinates z\ on this brane as z', = zijika^)- With these coordinates 
the warp factor in the metric (152j) is equal to one at the brane z = 02, and they are physical 
coordinates for an observer on this brane. The corresponding lengths of the compact dimensions 
are given by L\ = L\l(ka<i). For the part in the VEV of the fermionic current induced by the 
hidden brane and measured by an observer on the visible brane, from Eqs. (|45|) and (|53p one 
finds 



where the expressions for the primed quantities in the right-hand side are obtained from the 
corresponding expressions above by the replacement Li — > L',. When the left brane tends to the 
AdS boundary, a\ — > 0, the fermionic current goes to a finite limiting value. 

4 Fermionic current in finite-length carbon nanotubes 

In this section, we apply general formulas given above, for the investigation of the fermionic 
current in finite-length carbon nanotubes. Because of their unique electronic properties and 
potential applications in nanotechnology, carbon nanotubes have attracted enormous attention 
(see, for instance, Ref. [17J). A single- walled cylindrical nanotube can be thought of, as a 
graphene sheet rolled in a cylindrical form, with diameters ranging from 1 nm to 5 nm. Multi- 
walled carbon nanotubes are made of coaxial nanotube cylinders. The diameter of the outer 
tubes for multi-walled nanotubes can be of the order of 500 nm. The electronic coupling between 
separate layers in these nanotubes is weak, and the results given below can be applied for both 
single- walled and multi- walled carbon nanotubes. The low-energy excitations of the electronic 
subsystem in a graphene sheet are well described by the (2+l)-dimensional Dirac model con- 
taining a pair of spinors tpA and ips (f° r a review see Ref. [3]). The latter correspond to the 





71=1 



(54) 
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two different triangular sublattices of the honeycomb lattice of graphene. Recently, a number of 
predictions of the Dirac model have been confirmed experimentally with a high precision. For 
carbon nanotubes the background space for the (2+l)-dimensional Dirac model has the topology 
R 1 x S . The field equation for the separate spinors reads 

(m^V A> + il l Di ~ m)i>j = 0, (55) 

where J = A,B and vp ~ 10 8 cm/s is the Fermi velocity of electrons. The mass (gap) term in 
Eq. (|55p can be generated by a number of mechanisms (see, for example, Ref. }18j). 

The electrical properties of carbon nanotubes crucially depend on diameter and chirality. 
Depending on the chirality, the nanotubes are either metallic or semiconducting. The chirality 
also determines the periodicity condition along the compact dimension for the fields tpj. In 
metallic nanotubes, the periodic boundary condition is realized with the value cti = for the 
phase in Eq. ([3]). For semiconducting nanotubes, depending on the chiral vector, there are two 
classes of inequivalent periodicity conditions, corresponding to ai = ±1/3. These phases have 
opposite signs for the sublattices A and B. 

Although carbon nanotubes are typically longer than a micrometer, recently a number of 
techniques have been developed for the synthesis of ultra-short carbon nanotubes with lengths 
ranging between 20 and 80 nm (see, for instance, Ref. [19J ) . This type of nanotubes are es- 
pecially useful in biomedical applications. In the long-wavelength description of the electronic 
subsystem, the Dirac fields ipj live on the cylinder surface and for finite-length nanotubes ad- 
ditional boundary conditions should be imposed at the edges. As the electrons are confined 
between two edges of the nanotube, it is natural to impose the bag boundary conditions at the 
edges. As we mentioned before, the latter results in zero fermion flux through the boundaries. 
The effect of these boundary conditions on the fermionic current can be investigated by using 
the expressions from the previous section. For nanotubes the corresponding values for the di- 
mensions are as follows: D = 2, p = 0, q = 1. Another type of boundary induced effect, i.e. the 
Casimir effect between two graphene sheets or for a graphene sheet interacting with a metal, 
has been investigated in Refs. [20], using either the hydrodynamic model or the Dirac model 
for quasiparticles in graphene (for a comparison of the results obtained by these two approaches 
see Ref. [H]). 

Before considering finite- length effects, we recall the result for the fermionic current in an 
infinite-length nanotube. The corresponding expression is obtained from Eq. (|28p with a\ = a 
and L2 = L. Summing the contributions from the two sublattices and taking into account that, 
for them the phases a have opposite signs, for the total fermionic current we find [12] 

0' 2 >M = ^E cos (2^«)^(2vrncl>/$ )ii^, (56) 

n=l 

where a = 0, 1/3 for metallic and semiconducting nanotubes, respectively. In Eq. (|56p . we have 
expressed the component of the vector potential along the compact dimension in terms of the 
magnetic flux passing through the cross section of the nanotube: $ = A2L. Note that in Eq. 
(|56p and in the formulas below, we give the fermionic current for a given spin component. The 
total current is obtained multiplying by the number of spin components, which is 2 for graphene. 
From Eq. (|56p it follows that, in the absence of the magnetic flux, the total fermionic current 
vanishes, due to the cancellation of contributions from the two sublattices. 

Next, we consider a semi- infinite nanotube which corresponds to the single plate geometry 
discussed in the previous section. Summing the contributions from two different sublattices, 
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from the general formula (|33|) one gets: 



( ^(D 



(on) 



mvF 



E E 



with the notation 



777 



(±) 



j =+,-n= 



fc£ ±)2 + m 2 fc(±> 



-2xz 



CO 



dx 



(j)2 m + x' 



(57) 



2ir(n + <$>/<S> ±a)/L. 



(58) 



In Eq. (|57p . z = corresponds to the edge of the nanotube. An alternative expression is 
obtained by using Eq. ([65]) : 



/ 7 -2x(D 
V /(cn) 



mvp 

7T 2 



/•c 

■ V sin(27r($/$ +ia)) / 



x [cosh(L\/x 2 + m ) - cos(2?t($/$o + jot))] 1 

Z" 1 , 777 cos(2zxy) — xy sin(2zxy) 
x dy 
Jo 



(m 2 + x 2 y 2 )\/l — y^ 



(59) 



This expression is more convenient for the evaluation of the current near the edge and for 
large values of the nanotube diameter. The asymptotic expressions in various limiting cases 
are directly obtained from those for general D discussed in the previous section. The electric 
current corresponding to the VEV of the fermionic current discussed in this section is of the 
order \e\vpjL. The persistent currents in normal metal rings with this order of magnitude have 
been recently detected in Ref. [22] . 

In Fig. [3] we display the dependence of -^(j 2 )^) on ^ ne ma g ne tic flux, in units of the flux 
quantum, for mz = 0.5. The numbers near the curves correspond to the values of mL. The 
left/right panel correspond to metallic/semiconducting nanotubes. 




0.4 0.6 



Figure 3: Edge induced part in the fermionic current for semi-infinite metallic (left panel) and 
semiconducting (right panel) carbon nanotubes as a function of the magnetic flux. The graphs 
are plotted for mz = 0.5 and for different values of the parameter mL (numbers near the curves). 
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In the case of a nanotube of finite length, a, the general expression 

„-2> 



takes the form 



0' )(cn) - 0' )( cn ) ^ 

ri=— oo ^'=-1-,— 

1 + ^e a:c cosh [(a - 2z)x] 



U) 



dx- 



x±m 2ax i 1 



x 2 — ml J ^ 2 



(60) 



In the absence of a magnetic flux, the sublattices give opposite contributions to the fermionic 
current, and the latter vanishes. For a massless field one finds 



<Acn) = 0- 2 >[V^T £ £ ^E(- 1 )^o(2 S a|4 i) D- 



(61) 



n=— oo j=-f ,- 



8=1 



As already mentioned, the fermionic current is a periodic function of the magnetic flux, with 
period equal to the flux quantum $o- 

In Fig. U] we have plotted the boundary induced part, (j 2 )^) = 0' 2 )(cn) ~~ 2 )(cn)' f° r a 
massless field in metallic (left panel) and semiconducting (right panel) finite- length nanotubes, 
as a function of the magnetic flux. Numbers near the curves correspond to the values of the 
parameter a/L. 



a o.o 



V 

-J 




A 0.0 



V 



0.4 0.6 
<S>/<S> 




0.4 0.6 
<S>/<t> 



Figure 4: Boundary induced part in the VEV of the fermionic current for a massless field in 
finite-length metallic (left plot) and semiconducting (right plot) nanotubes, as a function of the 
magnetic flux. Numbers near the curves correspond to the values of the parameter a/L. 

The dependence of the current density on the length of the nanotube is displayed in Fig. 
[51 in the model with a massless fermionic field. The numbers near the curves correspond to 
the length of the compact dimension in units of a fixed length scale do, and for the magnetic 
flux we have taken the value $ = 0.8$o- The left and right panels correspond to metallic and 
semiconducting nanotubes. As it is seen, for long nanotubes the edge induced effects are small 
and the current tends the corresponding value for an infinite length nanotube. 

5 Conclusion 

We have investigated combined effects of compact spatial dimensions and boundaries on the 
VEV of the fermionic current. The geometry of boundaries is given by two parallel plates on 
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Figure 5: VEV of the fermionic current for a massless field in metallic (left plot) and semicon- 
ducting (right plot) nanotubes, as a function of the tube length. The graphs are plotted for 
<J> = 0.8^0 an d the numbers near the curves correspond to the values of the parameter L/clq. 

which the Dirac fermion field obeys MIT bag boundary conditions. Imposing these conditions 
in the region between the plates leads to the discrete eigenvalues for the component of the 
momentum perpendicular to the plates. For the case of a massive field, they are solutions of 
the transcendental equation (I12p . Along the compact dimensions, the periodicity conditions ([3]) 
are imposed with constant phases a\. In addition, we assume the presence of a constant gauge 
field which leads to the Aharonov-Bohm-like effect on the VEV of the fermionic current. The 
application of the summation formula (|15p allowed us to decompose the VEV as the sum of 
topological, single plate and interference parts. In this way the renormalization is reduced to 
the one for the boundary-free geometry. 

Similar to the case of the geometry without boundaries, the VEV of the charge density van- 
ishes. As for what concerns the spatial components of the current density, only the components 
along compact dimensions are non- vanishing. The VEV of the corresponding fermionic current 
depends on the phases in the periodicity conditions and on the components of the vector po- 
tential in the combination (|23p . The VEV is a periodic odd function of this parameter with 
the period equal to 1. It vanishes for special cases ai = and ai = 1/2. This is the case for 
untwisted and twisted fermion fields in the absence of the gauge field. 

Firstly, we have considered the geometry of a single plate. The boundary induced part in 
the VEV for the component of the fermionic current along the l-th compact dimension is given 
by Eq. (j33f) . This part vanishes for a massless field. An alternative expression, Eq. (j65[) . is 
derived in Appendix [A] by using the summation formula (|62p . Unlike the case of the fermionic 
condensate and the VEV of the energy-momentum tensor, the VEV of the fermionic current 
is finite on the boundary. In both limits of small and large values of the compact dimension 
length, the boundary induced part vanishes. At large distances from the plate, the asymptotic 
behavior of the fermionic current essentially depends on the phases present in the periodicity 
conditions along the compact dimensions. 

In the region between two plates, the VEV of the fermionic current is presented in the form 
(|42p . where the first term on the right-hand side is the corresponding VEV in the boundary- 
free space with topology RP +1 x (S l ) q . The second term is induced by the boundaries. An 
alternative representation is given by Eq. (|47p with the interference term defined by Eq. (|48p . 
For a massless field the distribution of the boundary induced part is uniform and is presented in 
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two equivalent forms, Eqs. and At large separations between the plates the boundary 
induced part is exponentially small (see Eq. ([15])), 

In the massless case the Dirac equation is conformally invariant for all spacetime dimen- 
sions. By taking into account that the bag boundary condition also is conformally invariant, 
and using the expression given above, we can generate the VEV of the fermionic current in 
conformally-flat spacetimes with compact spatial dimensions in the presence of boundaries. The 
corresponding VEV is given by Eq. (|53f) . where the first term in the right-hand side is the VEV 
for a conformally-flat spacetime with trivial topology, in the absence of boundaries, and the sec- 
ond term is induced by the nontrivial topology and by the boundaries. Important special cases 
correspond to background de Sitter and anti-de Sitter spacetimes. As an example, we have con- 
sidered the fermionic current induced on the visible brane in higher-dimensional generalizations 
of Randall-Sundrum-type braneworlds with internal spaces. 

As an application of the general formulas, in Sect. [H we have considered the VEV of the 
fermionic current in finite-length carbon nanotubes, described within the framework of the 
Dirac-like theory for the electronic subsystem. Recently a number of techniques allowed one 
to synthesize ultra-short carbon nanotubes which are especially useful in medical applications. 
For this type of nanotubes, the tube length can be of the order of magnitude of the compact 
dimension length (this is especially the case for the outermost tubes of multi-walled nanotubes) 
and the finite length effects may be important. The corresponding expressions for the VEV of the 
fermionic current are obtained by summing the contributions of the two triangular sublattices 
of the graphene lattice, with opposite signs of the phases along the compact dimension. Due 
to the cancellation between these sublattices, the current vanishes in the absence of a magnetic 
flux through the cross section of the nanotube. A magnetic flux through the cross section of the 
nanotube will break the symmetry allowing the current to flow along the compact dimension. 
The current is a periodic function of the flux with the period of the flux quantum <E>o- The 
pure topological contribution is then given by Eq. ([56]) . where the values a = and a = 
1/3 correspond to metallic and to semiconducting nanotubes, respectively. For a semi-infinite 
nanotube, the part induced by the edge is presented in two equivalent forms, Eqs. ([57]) and (f59|) . 
For a finite- length nanotube the VEV of the current is given by Eqs. (I60p and (|6ip for massive 
and massless fields, respectively. For long nanotubes the interference effects between the edges 
are small and the single edge parts dominate. 
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A Alternative representation for a single plate induced part 

An alternative representation for the part in the VEV of the fermionic current, induced by a 
single plate at z = 0, can be obtained by applying to the series over ni in Eq. (127ft the summation 
formula 123] 




oo 



«i=— oo 




A=±l 



(62) 
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For the corresponding series in Eq. (|27|) one has g{u) = u and the first integral in the right-hand 
side of this formula vanishes. By integrating over k p with the help of the analog of Eq. (|32l) . 
the VEV is presented in the form 

nr a i>oo / „2ixz „—2ixz \ 

(^)W = -^mL^naO £ / dx ( ^— + -L— - ) 
27rl/g Jo \m — ix m + ixj 



nq-ieZ^-i 

(y 2 -x 2 - m 2 )(p-i)/2 

— i,/- r s — 9 /q ~ \ ) (63) 

x z +m 2. cosh{yLi) — cos(2irai; 



with the notation 



a q-l 



D 

2 



^ [27rK + ai)/^] 2 +^ 2 - (64) 



Next, introducing instead of y a new integration variable i = yy 2 — x 2 — J7in g _ 1 , we pass to 

polar coordinates in the plane (x,t). After some additional transformations, the expression for 
the VEV of the fermionic current reads 



q n,_ 1 ez«-i J0 
x [cosh(L^/x 2 + m 2 ) — cos(27ra;)] _1 



1 m cos(2zxy) — xy sm{2zxy) 
^ (m 2 +x 2 y 2 )(l-y 2 )( 1 -P)/ 2 ' ( ' 

Note that in a model with a single compact dimension one has m Ilq _ 1 = m and the summation 
in Eq. (|65p is absent. By numerical calculations we have checked the equivalence of Eqs. f)33[) 
and (1651). 
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